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ABSTRACT 

We consider stationary, axisymmetric hydrodynamic accretion flows in Kerr 
geometry. As a plausible means of efficiently separating a small population of 
nonthermal particles from the bulk accretion flows, we investigate the formation 
of standing dissipative shocks, i.e. shocks at which fraction of the energy, angular 
momentum and mass fluxes do not participate in the shock transition of the 
flow that accretes onto the compact object but are lost into collimated (jets) or 
uncoUimated (winds) outflows. 

The mass loss fraction (at a shock front) is found to vary over a wide range 
(0% — 95%) depending on flow's angular momentum and energy. On the other 
hand, the associated energy loss fraction appears to be relatively low (< 1%) for 
a flow onto a non-rotating black hole case, whereas the fraction could be an order 
of magnitude higher (< 10%) for a flow onto a rapidly-rotating black hole. By 
estimating the escape velocity of the outflowing particles with a mass-accretion 
rate relevant for typical active galactic nuclei, we flnd that nearly 10% of the 
accreting mass could escape to form an outflow in a disk around a non-rotating 
black hole, while as much as 50% of the matter may contribute to outflows in a 
disk around a rapidly-rotating black hole. In the context of disk-jet paradigm, 
our model suggests that shock-driven outflows from accretion can occur in regions 
not too far from a central engine. 

Our results imply that a shock front under some conditions could serve as a 
plausible site where (nonthermal) seed particles of the outflows (jets/winds) are 
efficiently decoupled from bulk accretion. 



Subject headings: accretion, accretion disks — black hole physics — hydrody- 
namics — shock waves — galaxies: jets 
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Introduction 



It has been established by now that a large body of astrophysical objects hosting super- 
massive black holes (e.g., quasars, g alactic black hole candidates or rn icroquasars) exhibit 
collimated, powerful jets/winds (e.g.. iBegelman et al.lll984l : lLividll999l . for review). In par- 
ticular, strong outflows generally occur in the radio-loud active galactic nuclei (AGNs). It 
is now widely accepted that the outflows observed from a large class of objects as jets or 
winds have their origin in accretion flo ws, which at the sarne time power the radiatio n emis- 
sion associated with these objects (e.g.. iBlandford fc Payndll982l : iFender et al.ll2004l ). Many 
astrophysical systems apparently manage to transfer the energy from inflowing accretion 
flows to out f lowin g jets/winds of a small population of nonthermal particles. For instance, 
Junor et al.l (119991 ) in radio observations of the nearby active galaxy M87 found a remarkably 
broad jet with strong co Uimation occurr i ng alr eady at ~ 30 — 100 Schwarzschild radii from a 
central engine. Recently, iKataoka et al.l (120071 ) discussed the disk-jet connection of the radio 
galaxy 3C 120 observed with Suzaku. Such a population of relativistic outflowing particles 
is tho ught to produce a su bsequent synchrotron emission that is observed in several sources 



e.g., iMirabel et al.l 1 19991 ). This issue clearly points to the importance of investigating a 



fundamental connection between the accreting flows and outflows in regions not too far from 
the central engines - presumably within ~ 100 Schwarzschild radii. Because the emission 
of radiation requires the dissipation of the kinetic energy of the accretion flow, it is not 
unreasonable to suggest that the two phenomena, i.e. the dissipation/emission of radiation 
and the presence of outflows, are related to the same generic process, which manifests itself 
with different guises at the diverse sites that these phenomena are observed. 

Given that an accretion flow may very well become turbulent or, at a very minimum hot 
due to adiabatic compression, it is not surprising to anticipate such a flow to emit radiation, 
with luminosity and spectrum depending strongly on the speciflc character of the accretion 
process (e.g. quasi-spherical or disk). However, the ubiquitous presence of outflows in the 
same objects presents a different problem altogether: An outflow requires that a fraction 
of the accreted matter be endowed with velocity higher than the escape velocity associated 
with the speciflc radius at which the outflow is launched. Given that the only free energy 
available to the accreting gas is that of the gravitational fleld, hydrodynamic dissipation of 
the accretion kinetic energy can never produce an outflow since this process simply converts 
the kinetic energy of a gravitationally bound flow into thermal, with the energy per unit 
mass (thermal plus kinetic) never higher than that of the local gravitational potential. 

The launch of an outflow such as those observed requires the presence of an "engine" , i.e. 
a mechanism that expends mechanical work (i.e. low entropy energy) to transfer a fraction 
of the available energy to an even smaller fraction of the available mass, thereby imparting to 
it speciflc energy greater than the local gravitational potential; it is then expected that in its 
further evolution the excess energy will be converted into directed motion and the fraction 
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of the mass will escape to infinity in a coUimated (jet) or uncollimated (wind) outflow. 

The well known models of outflows, usually involving the action of magnetic fields such 
as the magnetocentrigual iet models (e.g..lBlandford fc Paynelll982l : IContopoulos fc Lovelace 
1994 : iKonigl fc Kartjdll994j : IVlahakis et al.ll200Cll ) , are specific examples as to what may con- 
stitute such an engine. These models assume the presence of a thin Keplerian disk "threaded" 
by a poloidal magnetic field; the Keplerian rotation of the magnetic field line footpoints and 
the magnetic tension transfer energy and angular momentum to the disk plasma which can 
escape to infinity. These models are consistent in that they solve simultaneously for the 
poloidal field geometry and the fl ow velocity (under c ertain simplifying assumptions, i.e. 
the self-similarity of the solutions). iPelletier fc Pudrita (119921 ) further developed the general 
theory of non-self-similar solutions of hydromagnetic disk winds. In these models the trans- 
fer of excess energy to the escaping particles is made at the expense of the rotational energy 
of the matter in the disk and it is mediated by the magnetic field. 

While this type of model gained popularity, the issue of jet/outlfow formation took a 
different turn with t he introduction of ady ection-dominated accretion flow (ADAF) (e.g., 
Narayan fc Yil ll994J : iManmoto et al.l 119971 ). These radiatively inefficient accretion flows 
(RIAF) were found to have positive Bernoulli integral of the flow and could therefore fulfill 
the condition necessary for the launching of jet/wind outflows; as such they present poten- 
tially interesting sites for the origin of such outflows. The po s itivity of the Bernoulli integral 
has been discussed and analysed by iBlandford &: BegelmanI (1l999l ) who pointed out that it 
is due to the combination of energy transfer by the viscous torques from the inner to the 
outer sections of the flow (the gas of the flow becomes bounded at its inner edge) and the 
local dissipation of the flow's azimuthal kinetic energy which is not radiated away but stored 
in the fluid to increase its internal energy. The latter authors then argued that the excess 
energy can be carried away to infinity (along with some fraction of the accreting mass and 
angular momentum) to produce continuous outflows from all radii to infinity while leaving 
the remaining flow with negative Bernoulli constant to naturally accrete onto the compact 
object: advection-dominated inflow-outflow solution (ADIOS). In this case, while the neces- 
sary excess energy is transferred by the viscous torques from the flow's more highly bound 
inner section, the necessary separation of mass to components with positive and negative 
total energy is still left unspecified. 

An al together different model th at offers a simplified picture of such a separation was pre- 
sented by ISubramanian et al.l (119991 ) who proposed that in the tenuous, collisionless plasma 
of an ADAF particles (protons) could be accelerated via a second-order Fermi acceleration 
by the shear motions of the underlying quasi-Keplerian azimuthal flow. They then argued 
that if sufficiently large pressure is built in the accelerated particle proton population (the 
electrons generally lose energy on time scales short compared to their transit time through 
the system and cannot build an energy density that could be dynamically important) and 
for favorable geometries of the disk magnetic field (large scale poloidal loops that open up 
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above the disk) the relativistic particle population could naturally (through the action of 
the gravitational field) segregate itself from the non-relativistic one, carrying off to infinity 
only the accelerated {E > rUpC^) portion of the disk plasma. In this case the engine is a 
combination of the particle acceleration and the action of the gravitational field. 



Finally, a model along the same lines was proposed by IContopoulos fc Kazanaj (119951 ) 
who suggested that even in the case of a completely turbulent magnetic field, a separation of 
the relativistic and non-relativistic particle populations is possible through the production of 
relativistic neutrons in the collisions of the relativistic protons with the ambient plasma and 
the ensuing production of relativistic neutrons. The subsequent decay of neutrons back into 
protons produces then a proton fiuid in regions of space devoid of inertia whose energy-to- 
mass ratio (and hence its asymptotic Lorentz factor) depends only on the ratio R/cTn, where 
R is the size of the system and the neutron life time and can lead to highly relativistic 
fiows for black hole masses M > lO^M©. 

In the present note we follow a similar simplified view to study outfiows in objects pow- 
ered by accretion: We consider the presence of (2-dimensional) shocks as a m eans of dissi- 
pation of the accretion kinetic energy in a fashion similar to that considered by IChakrabarti 
(jl990h and collaborators. That is we consider the transonic accretion of matter with the 
proper angular momentum to produce a standing shock at a radius close to the horizon, 
which, subsequently accretes onto the black hole after passing through a downstream sonic 
point. Previous steady-state analysis found that it is possible to judiciously choose the spe- 
cific angular momentum of the fiow, so that the outer transonic one could be connected 
through a shock transition to an inner transonic one which, passing through an inner sonic 
point, accretes onto the black hole. The location of the shock in this situation is determined 
by finding a radial position at which the density and velocity of the two fiow sections were 
those demanded by the dissipative Rankine-Hugoniot conditions across a shock. So far, the 
formation of standing shocks in hydrodynamic ac cretion has been extensively studied by 
a number of authors for both invisc i d fiows fe.g.. IChakrabartil Il990l: ISponholz fc Molteni 



Chakrabart 



1996; 



Lu &: Yuan 



1998 



Fukumura &: Tsurutal 1200J) and viscous fiows 



Chakrabartilll990l : iLu et al .I I1999I : IChakrabarti fc Dasll2004h . considering either dissi- 



1994 
(e.g., 

pative or non-dissipative shock jump conditions. Recently, standing shocks in the presence 
of poloidal magnetic fields [i.e., magnetohydrodynamic (MHD) shocks] around a black hole 
was also studied for various parameter dependence (see Das & Chakrabarti 2007 for pseudo- 
Newtonian geometry; Takahashi et al. 2002, Fukumura et al. 2007 for Kerr geometry). In 
the context of the particle acceleration via the first-order Fermi mechanism across a shock 
front, the production of shock-accelerated relativistic protons were discussed in spherical 
accretion (jProtheroe &: Kazanad Il983l : iKazanas fc EUisonl Il986f). while oth e r aut hors have 
explored the relativistic outfiows in ADAF with shocks (iLe &: Beckerl |2004| . |2005| ). Similar 
attempts have been made to make physical connections between the shocked-accretion and 
outfiows. For instance, mass outfiow rate were estimated from adiabatic shocked-fiow region 
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in Newtonian gravity (e.g. 
a similar approach to study the shock-generated outflows with little energy dissipation in 
pseudo-Newtonian geometry. Independently, from general relativistic MHD simulations, the 
formation of jets (magnetically-driven and gas-pressure dr iven iets) is found in the high pres 



sure regions due to the shock/adiabatic compression (e.g.. lKoide et al.lll999l : lNishikawa et al. 



20051 ). They concluded that the jets are mainly produced by the gas-pressure gradient which 
is greatly enhanced by the shock front at around r ~ 6 gravitational radii (note that this 
feature was not seen in the Newtonian calculations). These studies also suggest an essential 
connection between the shocked accretion flows and the jets; i.e., the shock front may serve 
base of the outflows. 

The novelty of our approach lies in considering the possibility of shock formation (i.e. 
obeying the general relativistic, dissipative Rankine-Hugoniot conditions at a shock front) 
in which part of the mass, angular momentum and energy fluxes escape in the z-direction 
and do not participate in the shock transition. We then examine the energy per unit mass 
of the escaping matter which we compare to the escape velocity at the shock radius; if 
it is greater than the latter we conclude that this scenario can produce an outflow with its 
outflow rate m and luminosity that are calculable and can be compared to those of the entire 
accretion to obtain a measure of the efficiency of our "engine" in producing outflows. Some 
population of energetic nonthermal particles, produced via a shock acceleration, may then 
be well separated from the equatorial accretion flows (which consists primarily of thermal 
particles) . 

More specifically, since we are interested in the formation of outfiows through shocks in 
the inner disk region relatively close to a presumably rotating black hole at a center (say, 
r < 30 gravitational radii), it is important to include the strong gravity and frame- dragging 
effects described by general relativity. To the best of our knowledge, no relevant work in 
the literature incorporates such mass and energy loss in the shock jump conditions. In the 
framework of our model presented here, mass loss is coupled to energy loss via the shock jump 
conditions, and therefore it must be considered simultaneously. It is this point that motivates 
us to explore, for the first time, the formation of outfiowing particles as a consequence of 
the formation of dissipative standing shocks in accretion in a fully relati vistic treatment. 
The f ormalism of our current model is partly based on the previous works (lYang fc Kafatos 
19951 : ILu fc YuanI 1 19981 : iFukumura fc Tsurutal I2004J ). Our main objective in this study is 
therefore to explore in details a possibility that the formation of outfiows (jets/winds) can 
occur at a dissipative shock front in transonic accreting fiows: i.e., a connection between 
shocked-accreting gas and the outfiowing particles. 

The structure of this paper is as follows. In §2 we review and explain our simplified model 
that simultaneously considers both shocks and outfiows with appropriate jump conditions. 
Main parameter dependence of the shock-outfiow solutions is explored in §3, where we show 
the nature of the shock-outfiow solutions and the corresponding global accretion solutions. 
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Our primary goal in these analysis is to examine the coupling between the shock-outflow 
solutions. In §4 we discuss our results and make some observation implications. Brief 
summary and concluding remarks are given there as well. 



2. Model Assumptions & Basic Equations 

In black hole accretion, accreting gas must be transonic. After passing through a first 
sonic radius, the gas is slowed down, and a shock may develop. For causality, the shocked 
gas must become supersonic again before crossing the event horizon. Below, we will explain 
the details of our simplified model. 



2.1. Accreting Flows around a Black Hole 

We consider a steady-state, axisymmetric accreting flows in Kerr geometry. The space- 
time metric is expressed by the Boyer-Lindquist coordinates as 

ds^ = - ( 1 — j dt^ dtd(p 

Zj Za 

where A = — 2mr + a^, E = + cos^ 9, A = (r^ + a^Y — a^A sin^ 9. m and a are 
mass and specific angular momentum (or Kerr parameter) of a black hole. Following the 
standard geometrized units, we have taken G = c = 1 in equation ([T]). Thus, the length 
(distance) r and a are measured in units of m. Since we are interested in the equatorial 
flows, we set = 7r/2 throughout this paper. The black hole horizon is then expressed by 
th = + yjm? — a^. The self-gravity of the flow is ignored, and we do not include the effects 
of magnetic fields for simplicity. 

We follow the earlier works on relativistic shock formation as follows; accretion time 
is assumed to be shorter than that of energy diffusion, thus treating the flows as adiabatic 
except at a shock front where a fraction of fluid energy, angular momentum, and mass are 
dissipated. To prescribe thermodynamic quantities we adopt a polytropic form as 

P = Kpl^^'^ , (2) 

where P and po are the thermal pressure and rest-mass density of the flow, which are locally 
measured in the fluid frame. Here, N denotes the polytropic index. The entropy of the fluid 
is characterized by K which is related to the entropy S hy S = log K where c„ denotes 
a specific volume heat of the flow. Because of stationarity {dt = 0) and axial symmetry 



{del) = 0), there exist two conserved quantities along a fluid stream line; namely, specific 
energy E and axial angular momentum L of the fluid, defined by 



E = -fiut , (3) 
L = fiu<p , (4) 

where /x = (P + p)/po is the relativistic enthalpy of the fluid, and p = Po + NP is the net 
baryon mass-energy density (including internal energy NP). The number density of the 
constituent baryon n is given by po = nrrip where rrip is the baryon mass in the flow. We 
assume that energy E and angular momentum L are both conserved along the flow except 
at a shock location (r = Vsh) where they are partially dissipated (will be explained in detail 
in the next section). 

From the four- velocity normalization of the flow {uaU°' = —1), we get 

l + M,M" + (M*)Vefr(r,A) = , (5) 

where the effective potential has been introduced by V^fj(r, A) = (7** — 2Xg'^'^ + X^g'^'^ with 
g^p i^ging the inverse metric components for g^p in equation ([T]). Here, the specific angular 
momentum of the flow A is defined by 

A-^ = -^. (6) 

E ut 

which is assumed to be conserved along the whole flow in our model. From equation ([5]) we 
obtain 



Ut{r,X) 



11 r H 1/2 

1 + UrU"^ 



-Kfr(r, A) 



From the definition of enthalpy and the polytropic relation given by equation (j2]) we 
can rewrite the enthalpy as 

p=l + {N+ l)Kpr • (8) 
Local adiabatic sound speed is defined as 



Combining the equations ([2]), ([8]) and ([9]) we can express p in terms of the sound speed Cg as 
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Accordingly, the energy of the flow E in equation ([3]) can be explicitly rewritten as a function 
of c. 



E 



1 + UrU''' 



cflf 



1/2 



/ (1 - Ncl) 



(11) 



which is conserved along the flow except at a shock front. That is, across a shock front, 
the energy E and the angular momentum L will both decrease in such a way that the ratio 
X = L/E is continuous across the shock. Using equation (jSj) and fllOl) the baryon rest-mass 
density can be rewritten as 



Po 



N 



[1 + 1/N){l-Nc^^)\ ' 
We define the mass-accretion rate M as 

M = -ATxrHvTpa , 
where H represents the vertical scale-height of the flow defined as 



H = 



(12) 



(13) 



(14) 



K 



from the conventional hydrostatic equilibrium assumption. Here, f^i<-(r) = m^/^/(r^/^ + 
am^/^) is the Keplerian angular velocity. Note that for accretion we have m** < 0. As 
often assumed, we only consider a constant mass-accretion rate i n this paper, although 



a var iable accretion rate has been discussed in the literature (e.g., iBlandford fc Begelman 



19991 ). Eliminating po from equations (fT2|) and ([T3l) we rewrite M as 



M 



c 



3.N+1 



{l + l/N){l-NclY 



(15) 



Similarly to E and L, the mass-accretion rate M is also conserved along the flow except at 
a shock location. After defining all the physical quantities necessary to solve for black hole 
accretion, we will describe below the transonic properties of the physical accretion solutions. 



2.2. Regularity Conditions 

For accretion to continue on to the event horizon, it is required that the accreting flows 
become supersonic at a sonic radius. After taking the derivatives of the equations ffTTl) and 
f|T5|) with respect to r (note that E and M are both constants), dcg/dr can be eliminated. 
Finally we obtain 
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where 



V = 2Ncl + Uru' {N {3cl - 2) - 1} , (17) 

^ dr 
- [a {ANcl + r=^/' {2Ncl (5 + 3^7„,) + v}] Kfr , (18) 



and 



V ^ mg,rcl + C{uY-^. (19) 
C = t{N{cI-2)-1] . (20) 

In equation ( |T6l) . the fact that the velocity gradient du^ /dr is finite at r = Tc requires 
that P(r = Tc) = and A/'(r = Tc) = simultaneously (i.e., regularity conditions), which 
will allow us to find the critical radius for a given parameter set. A physically valid 
accretion solution, therefore, must pass through a critical point at r = Tc before reaching 
the horizon whether or not the shock formation is possible. In the presence of a shock, a 
global shock-included accretion solution must go through a critical point on both sides of the 
shock location (i.e., before and after the shock forms). This fact requires multiple critical 
points. It is widely known that multiple critical points (up to three at most) can exist in 
general for a certain fiow parameter space. Much work has been done on examining the 
topological behavi ors (i.e., saddle, node, center and spiral points) of the critical points (e.g.. 



Chakrabartil Il990l ) . thus we will not repeat this. Because our main goal of this paper is to 
explore a possibility of outfiows that are coupled to the shocks in the global accretion fiows, 
we will only investigate the upstream fiows passing through the outer critical points r°"* 
while the downstream fiows passing through the inner critical points (the middle one is 
known to be unphysical). 

So far, for a specified fiow parameter set, one can obtain accretion solutions. Next, let 
us impose the shock conditions that connects one solution (i.e. upstream fiow) to another 
(i.e. downstream fiow), taking into account energy, angular momentum, and mass loss at a 
shock location. 



2.3. Shock Formation with Energy and Mass Loss 

Accreting fiows around a black hole are generally subject to a number of "invisible" 
obstacles that can decelerate the fiow: (1) centrifugal barrier due to the fiuid's angular 
momentum, (2) gas pressure- gradient, (3) radiation pressure-gradient, and (4) magnetic 
forces (i.e., pressure-gradient and/or tension force). Although our model is purely adiabatic 
and hydrodynamic [thus (3) and (4) are absent], fiows are still under the infiuence of the 
deceleration mechanisms (1) and (2). 
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Downstream 
Black Hole 



OutHows: 
Energy AE 

Angular Momentum AL 
Mass Outflow Rate AM 




Outflows 



Fig. 1. — Schematic picture of mass outflows originating from a standing shock front in 
hydrodynamic accretion. indicates azimuthal motion of the flow. At the shock, a fraction 
of the accreting gas is lost carrying energy AE, angular momentum AL and mass outflow 
rate AM. 



Following the previous works (see lYang fc Kafatoslll995l : lLu fc Yuanlll998l ). let us assume 
jump conditions that allow energy, angular momentum and mass loss at a standing shock 
front. Figure [1] illustrates a schematic description of our model. From equation ffTTl) we have 



El 




/ (1 - NCs2) 



(21) 
(22) 



where the subscripts "1" and "2" denote the quantities for upstream and downstream flows 
evaluated at a shock location (r = Vsh), respectively. We require Ei > E2 and define the 
energy dissipation and its fraction as 

AE = Ei- E2 and fE = ^, (23) 

where < Z^; < 1. The associated angular momentum carried by the outflows is then given 
by AL = Li - L2 = XAE. 

In the baryon mass conservation across the shock front, we also consider some mass loss 
(equivalently the loss of mass-accretion rate) associated with the outflows that are blown 
away as winds/jets. From equation f[T^ we have 



Ml = -47rr,h (d/2 + am'/^) u[ 



c 



2N+1 
si 



(1 + 1/iV) (1 - Ncl,) 



N 



Ki 



N ' 



(24) 



- 11 - 



M2 = -47rr3h (d/2 



am 



1/2 



Un 



C 



2N+1 
s2 



^2 



(25) 



.(1 + 1/7V)(1 

Clearly, the mass-accretion rate depends on the measurement of entropy K which must 
increase across the shock because of the heat generated (second law of thermodynamics). 
However, in the presence of both energy dissipation and mass loss at the shock front, a 
fraction of total energy (including mass and thermal energies) can be released from the flow 
surface, quickly reducing the rise of K. This hypothesis is justifiable when cooling processes 
are very efficient, although it is beyond the scope of this work to discuss the details of 
these mechanisms. Hence, for simplicity, we assume that the entropy will roughly remain 
unchanged at the shock front as a result from partial heat loss of the shocked flow. Thus, 
we set Ki = Ki = Kq and require Mi > M2. Similarly to energy dissipation, let us define 
the mass loss and its fraction as 



AM = Ml - M2 and fj^ = 



AM 

Ml 



(26) 



where < fj^ < 1. That is, fj^ = corresponds to no mass outfiows from shocks. 
The momentum fiux density T"'^ for an ideal fiuid is given by 

T'^f^ = (p + P)u''u'^ + Pg''^ . 
Therefore, the radial component T"^"^ is 



[I + l/N)Ur 

From the mass and momentum conservations in radial direction, we finally obtain 



M-lCs2 



--si 



'l + l/N)grrUl^ 



f^2Csl (1 - /m) 



;i + l/N)grrU. 



(27) 



(28) 



(29) 



which completes a series of our dissipative shock conditions. 



The strength of shocks is measured by the local compression ratio of the fiow, 712/711, 
which is expressed as 



n2 
ni 



^1 Csi| 



1 - /m) 



Cs2 



N 



(30) 



Hence, in the absence of mass loss [fg^j = 0), the compression ratio must be greater than unity 
when shocks occur. However, in the presence of mass loss (0 < fpj < 1), the compression 
ratio becomes a product of the velocity ratio and the mass loss fraction. Therefore, there 
can be a case where rarefaction or decompression with 7^2/7^l < 1 might take place even if 
the fiuid velocity (and sound speed) abruptly decreases across a very strong shock. To avoid 
further complications, we will focus our attention on compression shock waves only. 
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2.4. Dynamical Stability of Shocks 

In terms of the radial momentum balance across the shock front, some standing shocks 
can be dynamically unstable. That is, the shock may decay away (either radially inward or 
outward) as a result of a small (radial) perturbation of its position. In order to examine the 
stability of the obtained shock solutions, we perturb the radial momentum flux density T^"^ 
(equivalent to pressure) by invoking an infinitesimally small variation of the shock location 
(5rsh. If a shock front shifts back to its original location to retain the momentum equilibrium 
there, it is dynamically stable. The criteria is expressed as 



•rr 



dr 



dTl' 
dr 



5rsh = fi:(rsh)5rsh 



(31) 



sh 



where ^(rsh) is a function of r^h alone, thus possible to be numerically evaluated. In our 
calculations we take advant age of the known fact that ft:(rsh) < always guarantees the 
stabl e standing shocks (see 
2004h . 



Yang fc Kafatod Il995l : ILu fc YuanI 1 19981 : iFukumura fc Tsuruta 



3. Numerical Results 

Following the formalism outlined above we calculate the fractions of energy and mass 
losses for various model parameters. There are essentially three primary variables (E'l, A; r^h) 
that determine the solutions for a given geometry (a, 9) . In this paper, we restrict ourselves 
to the flows corotating with the black hole (aA > 0) at the equator {9 = 7r/2). 

The topology of the accreting solutions is normally classified as "x- type" or "g-type" 



depending on whether the shock-free solution is global or not (e.g., see iLu fc Yuan! 11998 
for definition). Although this is also one of the important aspects of the studies of stand- 
ing shocks, we do not consider the distinction here since it is not crucial to our current 
investigations. 

In our calculations we set = 3 and choose several representative values of the upstream 
flow energy: Ei = 1.003, 1.004, 1.005 for a/m = (Schwarzschild case) and 0.99 (Kerr case) 
to illustrate the frame- dragging effect. 



3.1. Dependence of Energy and Mass Loss 

We first present in Figures [2] and [3] the mass loss fraction as a function of shock 
location r^h for (a) Ei = 1.003, (b) 1.004, and (c) 1.005, with a/m = (Figure [2]) and 
0.99 (Figure [3]). Stable shocks are represented by solid curves while (dynamically) unstable 
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shocks are represented by dotted curves. Stable standing shocks, according to the criterion of 
equation ([ST]), can form in most cases in regions relatively close to the black hole {rsh/m < 
80). For a rotating black hole case, the shock location can be considerably closer to the 



hole (r^h/m 2 — 3), results similar to no mass loss cases (e.g., ISponholz &: Moltenilll994 



Lu fc Yuarull998l ). due to the fact that the horizon shifts more inward. There appears to be a 
smooth transition between the stable and unstable ones. Filled circles denote the maximum 
stable shock location, also corresponding to the weakest shock (i.e., smallest ^2/^1), while 
open circles denote the minimum stable shock location, also corresponding to the strongest 
shocks (i.e., largest n2/ni). In other words, the stronger dissipative shocks can develop 
at smaller radii, in agreement with the prev ious result in the absence of mass loss (e.g., 
Lu fc Yuanlll998l : iFukumura fc Tsurutal 12004 . in which /a^ = 0). Unstable shocks start to 
develop from where stable shocks become strongest. At both ends of the solution the curves 
are restricted by the dissipative shock conditions [equations (12T]) . (!22|) . (12^ . (125!) . fl29|) ] 
and the transonic properties [equation 0161) ]. It is seen that smaller angular momentum is 
required for shock (and accretion) to take place when the flu i d energy is larger, also consistent 



with previous studies of sta nding shocks (j Chakr abartill 1 99 6t iLu et al.lll997l : iLu fc Yuarull998 



Fukumura fc Tsurutal 20o3 ). Only stable shocks alone are allowed when angular momentum 
becomes sufficiently large (see Fig. ^ and Fig. [3]). To simplify our discussion, we will 
focus primarily on stable shocks alone from this point. As the angular momentum of the 
flow A increases, the shock location tends to (but not always) shift radially outward for a 
given energy Ei as expected because the centrifugal force correspondingly increases at a 
given radius. Hence, accreting flow subject to more outward force must decelerate at larger 
distance, forcing the shock location to shift outward, as seen in our results. This trend 
appears to be more clear for the a/m = 0.99 case. It is also noted that in these flgures larger 
angular momentum can allow for a larger mass outflow fraction fj^.^. 

In our model, mass loss fraction fj^^ can vary over a wide range (0% < f\,^ < 95%) de- 
pending on the shock location and angular momentum. That is, mass outflows are not sup- 
pressed by relativistic effects and energy dissipation (for comparion with psuredo-Newtonian 
case with no energy dissipation, see, e.g. Das &: Chakrabartil 1999 ). For a flxed angular mo- 
mentum, on the other hand, a higher value of f^j is expected from stronger shocks occurring 
in the inner regions, and this seems to be the case more in rotating black hole cases (see 
Fig. [3]). We shall explain the shaded regions in these flgures in the Discussion section. 

Figures m (a/m = 0) and [5] (a/m = 0.99) display the same mass loss fraction /j,;^ as 
a function of the energy loss fraction Z^;, corresponding to the solutions in Figure [2] and [3l 
respectively. At a flrst glance, there is an explicit positive correlation between fE and fj^j. As 
the angular momentum increases, maximum value of is clearly reduced. In these flgures 
the upper left portion corresponds to the outflow solutions with larger /^y ^^^d smaller (i.e., 
less energetic outflows) while the lower right portion represents the solutions with smaller 
fjyj and larger (i.e., more energetic outflows). Therefore, the results above suggest that 
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(a)E,= 1.003 (b ) E,= 1.004 (c)E,= 1.005 




10 15 20 30 50 10 15 20 30 50 10 15 20 30 50 

rsh/m hh/m hh/m 



Fig. 2. — vs. for various flow energies (a) Ei = 1.003, (b) 1.004, and (c) 1.005 with 
a/m = 0. Stable (unstable) shocks are denoted by solid (dotted) curves. In each branch, 
filled circles denote the maximum stable shock location also corresponding to the weakest 
shock (smallest n2/ni), while open circles denote the minimum stable shock location also 
corresponding to the strongest shocks (largest n2/ni). A labeled value of angular momentum 
A is fixed for each curve. Unbound outflow solutions are indicated by shaded boxes. 



{a)Ei= 1.003 (b)E, = 1.004 (c)E, = 1.005 




rsh/m rsh/m r^lm 

Fig. 3. — Same as Fig. [2] but for a/m = 0.99 case. 

more energetic outflowing particles may be separated from a shock front when the upstream 
flow possesses smaller angular momentum regardless of black hole spin a. Mass loss fraction 
/^j can become as high as < 95% for both a/m = and a/m = 0.99 cases regardless of 
energy Ei, whereas energy loss fraction is only fE < 1% for a/m = case but < 10% 
for a/m = 0.99 case. By direct comparison of these flgures, we note that rotation of a black 
hole a is also effective such that dfE/ df^,j{a/m = 0.99) > dfE/dfj^{a/m = 0), i.e. an increase 
in mass loss would allow more energy loss for a/m = 0.99 case. Therefore, more energetic 
outflows can be expected from a shock around a rotating black hole. 

The upstream flow energy Ei also affects the mass loss fraction. It is noted that the 
energy-dependence of fj^j is strong; higher energy Ei can lead to larger mass outflow fraction 
for a given angular momentum A and a shock location rsh. For instance, with A = 3.5 for 
a/m = case (in Fig. HI), < 20% when Ei = 1.003, 20 - 30% when Ei = 1.004 
and as large as ~ 40% when Ei = 1.005. With A = 2.16 for a/m = 0.99 (in Fig. [5]), 
/^y ~ 20 - 85% when Ei = 1.003, 40 - 85% when E^ = 1.004, and it is as high as 60 - 85% 
when El = 1.005. As a reference, we also examine the energy-dependence of the shock 
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Fig. 4. — vs. fE corresponding to the solutions in Fig. [2] for a/m = case. 
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Fig. 5. — Same as Fig. |4]but for a/m = 0.99 case. 



strength in Figure [6] where the (local) compression ratio ^2/^1 is plotted against Vgh for 
different energies: Ei = 1.003, 1.004 and 1.005 as used before. Angular momentum A is 
fixed in each case to see energy- dependence alone (although different values of A for different 
spin a must be chosen to obtain the solutions). We find that our dissipative shocks (coupled 
to mass loss) become stronger with decreasing energy, a behavior similar to the typical types 
of shock formation (e.g., see Lu et al. 1997 for adiabatic shocks; Lu & Yuan 1998, Fukumura 
& Tsuruta 2004 for isothermal shocks). 

To exclusively illustrate the black hole spin dependence a of the mass loss efficiency fj^^^ 
we fix all other parameters (i?i,A,rsh) except for a. Figure [7] shows f;^.^ against a and Je 
for A = 3.45 and Vsh/rn = 30. As seen in the earlier results, black hole rotation alone can 
clearly enhance the efficiency of mass outflows f^^ from ~ 3% (for a/m = 0) up to ~ 95% 
(for a/m = 0.35). On the other hand, the corresponding energy loss efficiency Je remains 
as low as ~ 0.02 — 0.1%. Note here that A would have to be properly adjusted in order to 
obtain the solutions for higher black hole spin a. 

We have chosen above some representative values for the flow energy for parametric 
purpose. Weakly viscous/inviscid accretion in general is a good model for some limited 
specific cases, like our Galactic center, for example. For such specific cases, the realistic 
choice of energy should be very small. From this perspective we examine to see whether low 
energy flows can still produce shock-driven outflows. Figure M shows mass loss efficiency fj(^ 
as a function of for a/m = 0. We set Ei = 1.000001 and A = 3.73. Mass outflows can 
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2 5 10 20 50 2 5 10 20 



Fig. 6. — n^jrix vs. rsh for (a) a/m = and (b) 0.99 cases. We choose E\ = 
1.003, 1.004, 1.005 from top to bottoms curves. Other notatfons are the same as before. 




02 0.4 

a/m 



Fig. 7. — Mass foss efficiency /^j as a function of black hole spin a and energy loss efficiency 
/e for a set of fixed parameters. We choose A = 3.45 and rgh/rn = 30. Solution curve is 
projected on each plane as shown. 

indeed be produced with fj^j ranging from ~ 1% up to ~ 65%. Both unstable and stable 
shocks are present as in the earlier cases but not continuously connected (no shock regions 
between the two). The range of shock location is much narrower in radius in this case, over 
which the mass loss efficiency can significantly change as mentioned above. We will discuss 
this more in the Discussion section. 

The major correlations we find among the primary parameters are summarized in Ta- 
ble [H Table [2] shows various correlations with shock strength. Compression ratio n2/ni is 
strongly correlated with (^'sh, /e, /a/)) stronger shocks are expected in regions very close to 
the central engine particularly around a rotating black hole. Strong shocks in principle are 
accompanied by high energy and mass loss fractions. 

To sum up, our results show that strongest stable shocks generally develop at the small- 
est radii (closer to the black hole), accompanied by the largest mass loss /^^^ (and the largest 
energy loss as well). The rotation of the black hole apparently amplifies the shock strength 



100 



no shocks 




6.5 7 7.5 8 8.5 9 9.5 10 



Tsh /m 

Fig. 8. — Mass loss efficiency fj(j vs. Vsh for a/m = 0. We set Ei = 1.000001 and A = 3.73. 
Notations are the same as in Figure [2J 

by more than a factor of two, also extending the outflowing site significantly inward (i.e., 
"^sh/^ ^2 — 3 when a/m = 0.99 while r^-^/m > 12 when a/m = 0). We will make some 
implications of the obtained shock-outflow solutions in the last section, §4. 



3.2. Global Accreting Flows 

Samples of the global run of the physical parameters of accretion flows that include 
shocks are given in Figures [9] [a/m = 0) and [10] {a/m = 0.99). In Tables [3] and H] we 
provide the specifics of the shocks associated respectively with the flows of Figures P and [TUl 
These plots have been made assuming that all flows accrete at 1% of the Eddington rate, 
i.e. m = M/MEdd = 0.01 and that the black hole mass is m = lO^M©, typical of Seyfert 
nuclei. Each panel shows (a) the radial velocity \u^\, (b) the angular velocity Q, (c) the 

Table 1. Correlation among primary quantities. 



Parameter 




fE 




A 


+ 




++ ^ 


El 






+ 



Note. — Strong positive (negative) cor- 
relation is denoted by H — h ( ) while rela- 
tively weak positive (negative) correlation is 
denoted by -|- (— ). No significant correlation 
is shown by <->. 
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electron scattering optical depth defined by r = narH, (d) flow temperature T [K], (e) the 
density po [g cm^^], and (f) the ratio of the vertical scale-height to the radius H/r. Vertical 
lines denote the positions of shocks connecting the upstream and downstream values of the 
corresponding quantities of each flow. Dotted curve in panel (b) shows the Keplerian angular 
velocity Qk- Note that for clarity purposes we only show three of the representative shock 
solutions obtained, although the shock can occur at any radius between the outermost and 
innermost solutions (with different values of /^;^ and Je)- 

We flnd that in both the non-rotating (Fig. [9]) and rotating (Fig. [TOl) cases the upstream 
flow density scales as po,i{r) ~ r"'^/^, as in the ADAF self-similar solution. On the other hand, 
as we explain below, the downstream flow density generally has a slightly steeper power-law 
slope, po,2{r) ~ r~'^/^ — unless the shock location is too close to the horizon. As seen 
in Figures [9^ and [TOk . this is primarily because of the decreasing (radial) downstream flow 
speed |m2('")I; whose feature is more obvious around a rotating black hole. After the shock 
transition, frame- dragging of the rotating black hole forces the accreting flow to accelerate 
more in the toroidal direction (i.e., corotate) rather than radial direction. Under adiabatic 
assumption, at the same time, the flow temperature continues to rise (Fig. [9li and [TOH) 
due to pdV compression of the flow. Since Cg oc T^/^, sound speed also increases. The 
resulting (thermal) pressure gradient {dP/dr < 0) prevents the incoming flow from speeding 
up and in fact decelerate the downstream flow (for a while) until the flow falls deep inside 
the gravitational potential well to radially speed up again. Consequently, the fluid motion 
is governed mainly by increasing toroidal velocity rather than radial one (compare Fig. [9^ 
with Fig. [TOk). This results in a larger radial velocity gradient d\u2\/dr > in Figure [TOh. 

In differentially-rotating flows, the Keplerian frequency increases with decreasing r even 
more rapidly than the sound speed, leading to a decreasing scale-height H [see equation (iMl) ]. 
The dependence of these two quantities, Qk and Cs, on the radius in the downstream flow 
could lead to geometrically slim/thick all the way to the horizon [see, e.g., H/r > 0.2 — 0.3 
in Figs. Ef and M]. 

Recalling po oc l/{rH\u^\) from equation f|T3|) . the above fact that both Iwgl and H 
decrease allows a very steep downstream density proflle Po,2('") after the shock transition 

Table 2. Correlation with shock strength ^2/^1 among primary quantities. 



Parameter 




fE 


f A/ 


712 /ni 


++ ++ 



Note. — Notations are the same 
as in Table [Tl 
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(and this is more so in a rotating black hole case because of larger drop in luj] as mentioned 
earlier). The downstream flows become slightly more opaque to electron scattering because 
of the shock compression, yet they continue to remain optically-thin (r < 1). 




Fig. 9. — Global shock-included solutions, (l)-(3), for a/m = case. See Table [3] for the 
parameters, (a) radial velocity \u^\, (b) angular velocity Q, (c) scattering optical depth r, 
(d) temperature T, (e) density po, and (f) aspect ratio H/r. Vertical lines denote shocks 
while filled circles indicate outer/inner critical radii. In (b) Keplerian profile Qx is denoted 
by a dotted curve. In (e) dotted gray lines show the slopes of —3/2 and —3. 



4. Discussion & Conclusions 

In this work we have examined the structure of accretion flows that includes shocks of 



more general character than those discussed so far in the literature (e.g.. IChakrabartilll990l ). 



Table 3. Sets of Parameters for the Global Solutions in Figure [9l 



Model 




n2/ni 


fE (%) 


/m (%) 


(1) 


15 


4.2 


0.319 


52.4 


(2) 


31 


3.8 


0.136 


47.3 


(3) 


54 


3.1 


0.00253 


43.4 



Note. — Note that a/m = 0, Ei = 1.003, A = 
3.6. 
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Fig. 10. — Same as Figure [9] but a/m = 0.99 case. See Table H] for the parameters. 

In particular we have examined whether it is possible to have flows with shocks at which 
part of the mass and/or energy fluxes are lost and do not participate in the shock transition 
(this maybe the case in multidimensional shocks or in shocks that involve acceleration of 
particles that escape from the shock region). Such generalized shocks obey jump conditions 
more general than those of Rankine-Hugoniot and it is not a priori certain if accretion flows 
can allow for such shocks. We have also examined the degree of mass and energy loss allowed 
that are at the same time consistent with the continuation of the downstream flow through 
another sonic transition onto the black hole. In this respect we have focused our study on 
shocks in which the energy per particle for those escaping the shock transition is greater than 
that of the local gravitational potential. In such a case one can argue that these particles 
can escape to infinity producing the jets/winds observed in many accretion-powered systems. 
We have found that there are indeed flows with global parameters (i.e. A and E) that allow 
for such shocks. On the other hand, we have also obtained solutions at which the energy per 

Table 4. Sets of Parameters for the Global Solutions in Figure [THl 



Model 




n2/ni 


fE (%) 


/m (%) 


(4) 


2.4 


10.2 


8.3 


89 


(5) 


10 


7.9 


2.6 


51.0 


(6) 


78 


2.73 


0.002 


17.6 



Note. — Note that a/m = 0.99, Ei = 
1.003, A = 2.16. 
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particle of the escaping matter is smaller than the local gravitational potential, indicating 
that in all likelihood these particles will stagnate and eventually join the rest of the flow 
onto the compact object. However, a deeper understanding of the precise mechanisms that 
can lead to the particles that do not participate in the shock transition, as we conjectured 
above, requires the knowledge of the detailed microphysics of the shock (i.e. the acceleration 
efficiency, the fraction of energy put into relativistic particles at the shock and their transport 
in and around the shock geometry); these are beyond the scope of the present paper. 

To get a crude estimate of the energetics of these outflowing particles, we compute a 
threshold energy of the escaping particles under a set of assumed parameters. The total 
energy flux carried by the escaping particles (of energy AE) is given by AE x ATcrshHnu^' 
[erg s~"^] where H = {Hi + H2)/2,n = (rii +n2)/2 and = {u{ + U2)/2 are evaluated at the 
shock location (r = Vsh)- The minimum kinetic energy necessary for the separated particles 
(of mass outflow rate AM) to escape the bulk flow is given by AMv'^/2 [erg s~^] where Vc, 
to be determined below, is the required critical velocity of the particles. By equating these, 
we solve for Vc to obtain 

~ 8nr,^Hnu'~mpC^AE/AM . (32) 

For the outflowing particles to be effectively decoupled from the bulk flow, kinetic energy 
associated with this critical velocity should exceed or be at least comparable to the corre- 
sponding gravitational binding energy at the shock location. This condition reads 

2"' ~ • 

or 



AE>E,, (34) 

where Ec = GmAM / {Anr'l^HnvJ'mpC^). Since all these parameters {H,n,u^) are functions 
of shock location r^h which depends on /^j, Ec is obtained once we specify fj^^. Now, we can 
review the shock solutions obtained earlier from this new perspective: In Figures [2H5] the 
allowed shock solutions satisfying AE > E^ which is relevant for outflows that may escape 
the shocked flow, are displayed by shaded regions, while the rest of the solutions correspond 
to the outflows presumably bound to the bulk flow. According to these estimates, an outflow 
that expels a certain amount of the preshock accreting matter is possible, over a rather broad 
range of radii in the vicinity of the central object; up to nearly 10% of matter may participate 
in producing seed particles for jets/winds within ~ 20 — 40 gravitational radii around a non- 
rotating black hole, while as much as 50% of preshock matter could contribute to forming 
jets/winds within ~ 2 — 30 gravitational radii around a rapidly-rotating black hole. More 
interestingly, our solutions also allow for dissipative shocks {Je 7^ 0) with negligible mass 
loss {fj(j ~ 0). Such solutions that involve the loss of inflnitesimal mass and flnite amount of 
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energy can lead to relativistic outflows from the corresponding shocks. Perhaps, outflows via 
this type of shocks may produce kinematically strong jets (with high Lorentz factor) rather 
than moderate velocity winds. Thus, shock-driven outflows may provide clues on the origin 
of the jets/winds from the inner accretion regions, for instance, in some active galaxies, e.g., 
M87 and 3C120. 

In this simple model, shocks in principle could drive axisymmetric outflowing matter. 
Near the base of the outflow, as illustrated in Figure [H its shape would be that of a hollow 
cone. However, as it propagates to longer distances we do not expect that this shape would 
be preserved all the way. Within the framework of our current scenario it can be speculated 
that the most of the outflows discussed here could be in a diffuse form, while some strong 
outflows with low fj^j and high Je might possess a relatively more collimated geometry 
(perhaps in the poloidal direction) due to its high Lorentz factor. In the presence of large- 
scale magnetic fields the dissipated plasma would stream along the field lines, producing 
collimated outflows. However, it is beyond the scope of our model to further discuss the 
exact geometry of the outflows. 

It should be noted that for a given A there exists the degeneracy of shock-outflow 
solutions: i.e., for a single value of angular momentum A there is a finite range of possible 
shock locations allowed (r^™ < r^h < f'^^) that corresponds respectively to outflows with 
f^"^ < /m < f^^^ where the indices "min" and "max" denote minimum and maximum 
values, respectively. In the context of our model, one can predict the shock location by 
the (observational) knowledge of (A, Ei, /^j), although the obtained parameter space can be 
altered by additional physical ingredients (e.g., viscosity, magnetic fields, for instance). 

Another issue to be addressed is the fact that in general critical points are not necessarily 
equivalent to sonic points depending on the flow geometry and equation of state used (e.g.. 



Dad 120071 ). That is, there could be a potential danger that a false shock could occur in 
subsonic region between a critical radius and an actual sonic radius, in which case the 
obtained shock would be unphysical. To ensure that a physically valid shock forms in 
supersonic regions, we have checked the validity of our shock solutions by computing a three- 
velocity component of the flow measured by a suitable local observer at the shock location. 
We calculated the (radial) flow velocity v " ^ measured by a l ocally stationary observer in the 
corotating reference frame (e.g.. lLulll986l : ILu fc Yuanlll998l ) 



1 + UrU^ 



(35) 



and compared this velocity to the local sound velocity Cg given by equation Q. All the shock 
solutions presented here are found to form in supersonic regions (i.e., \v'^\ > at r = rsh). 

We showed that low energy flows can still produce mass outflows with suitable angular 
momenta. Although continuous accreting solutions (i.e. shock-free solutions) are persistently 
present even for smaller energy Ei (as expected), we do not find shocked flow solutions (or 
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perhaps they are present but in much narrower parameter space). It can be speculated that 
maybe by prohibiting our jump condition in energy (energy dissipation at a shock front) we 
may obtain shock-d r iven o utflows for smaller flow energy Ei , which could be a case similar to 



Das &: Chakrabartil (119991 ) where they considered little energy loss and found mass outflows 
in pseudo-Newtonian geometry. However, powerful outflows should carry away a signiflcant 
amount of (kinetic) energy. Hence, outflow solutions should in principle be coupled to energy 
dissipation as treated here. 

We have explored a coupling between shock solutions in accretion and mass/energy 
losses (fractions) under a scenario that the shock-driven outflowing particles may participate 
in forming a base of jets/winds. For various flow parameters with a given black hole spin, we 
have shown, by steady-state, axisymmetric hydrodynamic calculations, that the dissipative 
shock front could be a plausible site where a fraction of the accreting matter can be decoupled 
as jets/winds from the bulk accretion flows. 
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